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Introducing fermionic loops contributions in Numerical Stochastic Perturbation Theory was mainly motivated 
by the proposal to compute 2-3 loops for renormalization constants (and improvement coefficients). This is feasible 
because the computational overhead of unquenching NSPT is by far lower than for non perturbative simulations. 
We report on first, preliminary results for the quark propagator (basically the third loop for the critical mass) 
and discuss the status of the computation of other quantities. 



INTRODUCTION 

Numerical Stochastic Perturbation Theory 
(NSPT) is a numerical method to perform Lat- 
tice Perturbation Theory which has proven to be 
very powerful in the quenched approximation of 
Lattice QCD. In recent years PQ the unquenched 
version was introduced, whose first (benchmark) 
results were encouraging The goal, which 

appears feasible, is to have NSPT competitive 
on any 2-3 loops computation of quantities like 
renormalization constants (and possibly improve- 
ment coefficients, where needed) for potentially 
any fermionic action. An obvious first benchmark 
one can go through is the quark propagator for 
Wilson action. In other words, one is comput- 
ing to 3 loops, for example matching to the 
RF scheme. In the following we only report on a 
byproduct, the critical mass for Nf = 2. 



1. The computational set-up 

For a comprehensive account of the formulae of 
unquenched NSPT we refer the reader to £Q. We 
briefly discuss instead the computational set-up. 
Results are obtained on a 32 4 lattice with Nf = 2 
Wilson action. The relevant configurations are 
generated on a APEmille crate (128 FPU's for 64 
GFlops of peak performance). Smaller lattices to 



keep finite size effects under control are demanded 
to a PC cluster. Stability with respect to single vs 
double precision has been tested. Gauge is fixed 
(where needed) to the Landau condition (Fourier 
acceleration has turned out to be very helpful). 
Given the replication of the fields at every per- 
turbative order, a single configuration on a 32 4 
lattice at 3 loops (i.e. f3~ 3 ) level amounts to 1.8 
GB. We are collecting a large amount of confi- 
gurations in order to set up a repository out of 
which any fermionic observable can be computed 
in perturbation theory to 3 loops. 

1.1. Status of unquenched computations 

As already stated, the configurations we have 
been generating up to now are at Nf = 2. Given 
the polynomial dependence of fermionic observa- 
ble on Nf in perturbation theory, it will be easy 
to fit it once configurations will be available also 
at other value of Nf. There is anyway a computa- 
tion already on its way just now, which is that of 
the third order of the residual mass counterterm 
in Lattice Heavy Quark Effective Theory. The 
goal is to extend the unquenched result for the 
b-quark mass [3] to three loops accuracy for the 
perturbative subtraction which is needed to keep 
renormalons effects under control. An Nf = 2 
computation will be enough for this purpouse. 
In recent years other methods to compute the b- 
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Figure 1. —So at order /3 _1 as a function of p 2 
(Nf = 2) in renormalized perturbation theory 
(critical mass counterterm as known from ana- 
lytical computations plugged in at one and two 
loops level). 

quark mass in completely non-pcrturbative ways 
have been proposed @|- The only unqucnched 
computations available at the moment are any- 
way still the ones relying on HQET and perturba- 
tive subtractions of the residual mass, so we still 
think that it is worthwhile to proceed. The resid- 
ual mass, as it was the case for the quenched case, 
is being computed out of the computation of the 
matching between lattice and potential scheme 
couplings, which in turns relies on computations 
of big Wilson loops. Quark bilinears are the ob- 
vious first quantities we are computing (again, at 
Nf = 2 level, at the moment) in the quark sec- 
tor of the theory. As a matter of fact also four 
fermions operators look feasible. 

2. The critical mass 

The quark propagator (which we denote by 
S(p ,m)) is our first benchmark computation 
among the quark bilinears. In NSPT approach 
one actually computes S and then proceeds to 
invert to get 

T 2 (p 2 ,m) = Sip^m)- 1 (1) 
= i — m — S(p 2 , to) 



which is the quantity which actually encodes all 
the information. S in turns goes like 

S(p 2 , to) = S + i jffii(p 2 , to) + toS 2 (p 2 , m) (2) 

where So is a pure lattice artifact (actually di- 
verging as a -1 ). At to = it is the well-known 
critical mass, the additive renormalization one 
has to subtract from the bare lagrangian mass 
in order to obtain the bare mass that gets multi- 
plicatively renormalized. It is the only quantity 
we are actually going to report on here (again, for 
Nf = 2). Still it is worthwhile to note the fea- 
tures all Si share with respect to our NSPT com- 
putation. As already stated, one first computes 
S 1 , then inverts to get T and then projects on the 
(Dirac space) identity and the 7 M in order to get 
the S^ Of course one takes advantage of the 
colour simmctry to improve the accuracy. The 
dependences on p 2 and to are completely known 
at one loop level, while those on p 2 and to = 
are known from continuum perturbation theory 
(notice that from these one can elaborate on p 2 
and to dependence at 2 and 3 loops). The to = 
results are the ones we are interested in here (they 
also define the point where one would like to com- 
pute in general the renormalization factors Z). 
We used two IR regulators in order to get the 
to = results. We both plugged small m « val- 
ues in and subtracted zero modes in the to = 
simulations. Results showed that we can control 
the procedure to a good accuracy. On the lattice 
at a finite value of the lattice spacing S ^ S(p 2 ). 
The dependence on the lattice invariants one can 
construct as powers of pa are actually the natural 
handles for the a — + limit. 

The first and second loop contribution to the crit- 
ical mass are analytically known 0. For Nf = 2 
they read 

- S = 2.6057 /3~ 1 + 4.293 /T 2 + . . . (3) 

Notice that once one knows these values, one ac- 
tually has to plug them in as counterterms in 
NSPT simulations (to keep their contributions 
into account at higher loops). In Fig. ^ we plot 
the component along the identity of T 2 (i.e. So, 
given the to = value) at order This is 

plotted against p 2 . The curve is not smooth since 
(as we have already pointed out) we know this is 
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Figure 2. — Eo at order (3~ 2 as a function of p 2 
(Nf = 2) in renormalized perturbation theory 
(same framework as in Fig. ^) . 
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Figure 3. — Eo at order /3 -3 as a function of p 2 
(Nf = 2). The result has not yet been extrapo- 
lated in the stochastic time-step r — * 0. 



not a function of p 2 alone: the value extrapolated 
at p = is the quantity one is interested in. No- 
tice that this is what one would call renormalized 
(with respect to the additive mass renormaliza- 
tion) perturbation theory: one has plugged the 
first two counterterms in, so Eo gets cancelled at 
one and two loops level. The confirmation of this 
at two loops level is clear in Fig. [5] (the same as 
Fig. [I] at I3~ 2 order). From Fig. one can in- 
spect that the result for the third order is within 
reach. We are not quoting a number at this pre- 
liminary stage, but the order of magnitude of the 
result is clear. Notice, by the way, that we plot 
the result from a very reduced set of configura- 
tions at a fixed time-step (in NSPT one performs 
simulations at a fixed value of the time-step used 
in integrating Langevin equation and then has to 
extrapolate to the continuum-time stochastic dy- 
namics). We add just a few comments to have a 
rough idea of the order of magnitude of this new 
contribution to the (perturbative) critical mass. 
If one compares the perturbative result for the 
critical mass at three loops level with the one 
determined with non-perturbative methods, one 
finds that results are still quite far away from each 
other. To have an idea of the quite poor conver- 
gence of the series in the bare coupling: at f3 = 5.5 
the contribution of order /3 -3 is some 50% that 



of order (3 2 (which is in turns 30% that of order 

Conclusions and acknowledgments 

By performing a three loops computation of the 
Wilson action critical mass at Nf = 2 we showed 
that three loops computations are indeed feasible 
in unquenched NSPT. The goal is now to compute 
at this order various renormalization constants. 
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ful discussions. 

REFERENCES 

1. F. Di Renzo and L. Scorzato, Nucl. Phys. B 
Proc. Suppl. 94 (2001) 567. 

2. F. Di Renzo, V. Miccio and L. Scorzato, Nucl. 
Phys. B Proc. Suppl. 119 (2003) 1003. 

3. V. Gimenez, L. Giusti, G. Martinelli and F. 
Rapuano, JHEP 03 (2000) 018. F. Di Renzo 
and L. Scorzato, JHEP 02 (2001) 020. 

4. R. Sommer, Nucl. Phys. B Proc. Suppl. 119 
(2003) 185. G. M. de Divitiis, M. Guagnelli, 
R. Petronzio, N. Tantalo and F. Palombi, 
|hep-lat/0 305018. 

5. E. Follana, H. Panagopoulos, Phys. Rev. D 
63 (2001) 17501; S. Caracciolo, A. Pelissetto, 
A. Rago, Phys. Rev. D 64 (2001) 94506. 



